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INTRODUCTION 


In  recent  years,  considerable  progress  has  been  made  in  the  development 
of  computer  programs  to  perform  dynamic,  large  deflection,  elasto-plastic 
analysis  of  submarine  structures  (1,  2,  3).  A  typical  problem  of  interest 
is  the  dynamic  elasto-plastic  analysis  of  submerged,  stiffened  cylinders. 
Stiffeners  on  such  cylinders  are  usually  thin  walled  open  cross  sections 
such  as  T,  Z,  L  or  flat  bar  sections. 

It  is  well  known  that  if  a  section  has  an  enforced  center  of  rotation, 
(i.e.  it  is  not  free  to  twist  about  its  shear  center,  as  a  stiffener  welded 
to  a  shell  is  not),  there  is  an  increase  in  the  warping  resistance  of  the 
section  to  torsion  (13,  22).  For  certain  cases,  it  has  been  shown  that  the 
stresses  resulting  from  resistance  to  warping  can  have  considerable  effect 
on  member  response  and  stability  (12,  13,  23).  In  addition,  experiments 
(7,  8,  19)  have  clearly  indicated  that  importance  of  stiffener  tripping 
(lateral-torsional  instability)  as  a  critical  failure  mode  for  ship  structures. 
Hence,  it  becomes  desirable  to  consider  restrained  warping  stresses  in  the 
development  of  an  elasto-plastic  stiffener  for  inclusion  in  the  analysis  of 
stiffened  cylinders. 

The  inclusion  of  warping  stresses  in  the  elastic  response  of  thin 
walled  open  cross  sections  has  been  studied  by  many  authors  in  great 
detail  (13,  14,  22).  The  plastic  response  of  the  section  under  combined 
stresses  is  defined  by  a  yield  surface.  Following  procedures  suggested  by 
Hodge  (16),  Morris  and  Fenves  (19)  have  determined  an  approximate  lower 
bound  to  the  exact  yield  surface  for  a  variety  of  cross  sectional  types. 

This  investigation  included  Saint  Venant  torsion  but  neglected  warping 
restraint.  Boulton  (5),  Dinno  (9,  10)  and  Gjelsvik  (13)  have  studied  the 
effect  of  warping  restraint  on  plastic  collapse  of  thin  walled  open  cross 
sections,  but  only  for  simplified  cases  in  which  one  or  more  of  the  possible 
stress  resultants  are  zero. 
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This  report  presents  the  formulation  of  an  elasto-plastic  thin  walled 
open  cross  section  when  the  effect  of  warping  torsion  is  more  important 
than  Saint  Venant  torsion.  Equations  are  developed  for  the  yield  surface 
of  a  Z-section  in  terms  of  four  stress  resultants,  namely  the  axial  force, 
bending  moments  in  two  planes,  and  the  warping  moment.  These  are  shown  to 
satisfy  both  upper  and  lower  bound  theorems  of  plasticity  and  therefore 
correspond  to  exact  yield  surfaces  within  the  assumptions  of  thin  walled 
sections . 

The  development  of  a  yield  surface  including  warping  restraint  is 
a  necessary  first  step  in  the  analysis  of  inelastic  stiffener  tripping. 

The  equations  are  currently  being  incorporated  into  the  computer  program 
EPSA  (2),  developed  at  Weidlinger  Associates  under  DNA/ONR  sponsorship. 

This  implementation  is  aimed  at  extending  the  modeling  capability  of  EPSA 
to  include  stiffener  tripping  and  improving  structural  response  predictions 
in  critical  regions  of  large  strain  and  deformation.  This  capability  will 
be  required  in  the  analysis  of  a  series  of  plastic  range  tests  on  both 
torpedo  (7)  and  submarine  scaled  model  (8)  structures. 
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KINEMATICS  OF  THIN  WALLED  OPEN  CROSS  SECTIONS 

The  assumptions  made  in  the  formulation  of  the  theory  of  thin 
walled  beams  are  that  cross  sections  deform  as  rigid  bodies,  and  that 
the  shear  strain  at  the  middle  surface  of  the  section  is  zero  (13,  22). 

Referring  to  Fig.  1,  the  axial  displacement  of  a  point  on  the 
middle  surface  of  the  bar,  w(s,  z) ,  may  be  written  as  (13), 

w(s,  z)  *  W(z)  -  U'(z)  x(s)  -  V '  (z)  y(s)  -  <j>'( z)  (u(s)  (1) 

where  n-s-z  are  coordinates  with  n  normal  to,  and  s  following  the  con¬ 
tour  of  the  middle  surface;  U(z),  V(z),  W(z)  represent  the  translation 
of  the  center  of  rotation  P  in  the  x-y-z  coordinate  system;  4> ( z)  is 
the  rotation  of  the  section  about  the  pole  P  ;  u>(s)  is  the  warping  func¬ 
tion  defined  as 

ui(s)  =  |  r(s)  ds  (2) 

o 

where  r(s)  is  the  perpendicular  distance  form  the  point  P  to  the  point 
in  question  on  the  middle  surface.  Thickness  or  secondary  warping  is 
neglected  in  Equation  2. 

The  first  three  terms  on  the  right  hand  side  of  Equation  1  satisfy 
Navier's  assumption  while  the  last  term  represents  warping. 

Stress  and  Strain  Components 

The  only  non-zero  stress  components  are  the  longitudinal  stress, 

0  ,  and  the  shearing  stress,  T  ,  shown  in  Fig.  2.  All  stress  and 

zz  zs 

strain  components  normal  to  the  wall  are  neglected.  In  addition,  all 
in-plane  strains  are  assumed  small.  From  Eq.  1,  the  longitudinal  strain 
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IS 


£  (s,z)  =  W’(z)  -  U"(z)  x(s)  -  V" (z)  y (s)  -  4>"(z)  lo(s)  (3) 

zz 

If  it  is  desired  Co  extend  this  formulation  to  account  for  geometric 
nonlinearities,  the  appropriate  nonlinear  strain-displacement  relations 
would  have  to  be  introduced.  This  has  been  done  in  (6). 

Stress  Resultants 

The  strain  energy  per  unit  length  of  the  bar,  l’s>  can  be  expressed 
as  (13), 

2U  *  NW'  +  M  U"-  M  V"  +  M  <$>"+  T  d>  ’  (4) 

s  y  x  u)  s 

In  Eq.  4,  the  stress  resultants  N,  M  ,  and  M  represent  the  axial  force 

y  ^ 

and  bending  moments  about  the  y-axis  and  x-axis,  respectively.  Saint 

Venant  torque  is  denoted  by  Tg,  and  is  the  warping  moment,  or 

bimoment,  and  is  defined  as 

M  *  -  a  JJ  dA  (5) 

U)  zz 

A 

These  stress  resultants  are  shown  in  Fig.  3. 
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DECOUPLING  OF  THE  MIXED  TORSION  PROBLEM 

The  last  two  terms  of  the  right  hand  side  of  Eq.  4  are  the  contribution 
of  the  warping  moment  and  Saint  Venant  torsion  to  the  strain  energy  of  the 
bar.  In  certain  applications,  effects  of  one  of  these  may  predominate.  In 
the  elastic  range,  the  relative  importance  of  Saint  Venant  torsion  to  the 
warping  torsion  in  straight  beams  is  given  by  the  magnitude  of  a  dimension¬ 
less  bar  parameter,  u^,  defined  in  (13)  as 

\  =  L/d  (6) 

in  which  L  is  the  span  length  of  the  bar  and  d  is  a  characteristic 
length  defined  by 

d2  =  El  /GJ  (7) 

IlXjJ 

where  E  and  G  denote  Young's  modulus  and  the  shear  modulus  and  I 

<j&i 

and  J  are  the  warping  constant  and  torsion  constant. 

It  has  been  shown,  (13  ,  23,  that  when  y^  is  small  ^approximately 

1)  warping  torsion  is  dominant  and  the  Saint  Venant  torsion  term  in  Eq.  4 
may  be  neglected.  If  y^  is  large  (approximately  20) ,  Saint  Venant  torsion 
controls  and  the  term  containing  the  warping  moment  in  Eq.  4  may  be 
dropped. 

In  the  plastic  range,  the  relative  importance  of  these  two  effects 
may  be  estimated  by  using  Eq.  6  with  d  assumed  equal  to  its  elastic 
value  and  L  taken  as  an  effective  length  corresponding  to  plastic 
deformation  patterns  of  the  bar. 

Case  Study 

The  preceding  discussion  on  approximations  to  the  mixed  torsion 
problem  is  applied  to  a  typical  problem  of  interest.  Consider  the 
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ring-stiffened  cylinder  in  Fig.  4.  A  cross  section  of  a  typical  stiffener 
is  also  shown.  The  value  for  the  warping  constant,  I  ,  for  the  2- 

LuUj 

stiffener  with  an  enforced  center  of  rotation  about  the  line  of  attach- 

6  6  6 

ment  between  the  stiffener  and  the  shell  is  .0176  in  (4.73  x  10  mm  ) . 

4  4 

The  value  of  the  torsional  constant,  J,  is  .00253  in  (1053  mm  ),  (14,  15). 

2 

Using  a  value  for  E/G  =  2.6  for  the  steel  structure,  d  given  by  Eq.  7 
is  18  in“  (11610  mm^). 

It  has  been  shown  experimentally  and  analytically  in  (4,  18)  that 

for  cylinders  with  relatively  stiff  ring-stiffeners,  plastic  deformational 

patterns  are  characterized  by  a  large  number  of  circumferential  waves 

(typically  greater  than  five).  A  conservative  estimate  of  an  upper  bound 

to  the  effective  length  of  the  ring-stiffener  may  therefore  be  taken  to 

be  that  corresponding  to  the  n  =  5  deformational  mode.  The  resulting 

value  of  should  also  be  an  upper  bound.  The  effective  length  is  then 

one-tenth  of  the  ring  circumference.  Using  a  value  of  8.9  in  (226  mm) 

for  the  length  of  the  stiffener,  a  value  of  *  2  is  obtained.  With  this 

value  of  as  an  upper  bound.  Saint  Venant  torsion  mav  be  neglected  in 

Eq.  4,  and  only  the  four  stress  resultants  contributing  to  the  axial 

stress,  0  ,  need  be  considered, 

zz 


PLASTIC  BEHAVIOR 


In  the  formulation  of  a  plastic  theory  of  thin  walled  bars,  the 
assumption  is  made  that  the  kinematics  describing  the  plastic  behavior 
of  the  bar  are  the  same  as  in  the  elastic  range,  (13)- 

The  plastic  response  of  a  bar  is  defined  by  a  yield  surface. 

The  assumptions  made  in  the  derivation  of  this  surface,  in  addition 
to  the  kinematic  ones  are  the  following: 

1)  the  material  is  elastic-perfectly  plastic, 

2)  the  von  Mises  yield  condition  applies, 

3)  shear  stresses  due  to  warping  are  ignored  as  in 
plastic  beam  theory,  (20}, 

4)  the  yield  surface  is  a  closed,  convex  surface  and  the 
incremental  plastic  deformation  is  directed  outward, 
normal  to  the  surface.  With  these  assumptions,  Drucker's 
hypotheses, (11) ,  are  satisfied  insuring  both  uniqueness 
and  stability  of  plastic  response. 

The  validity  of  any  proposed  yield  surface  may  be  assessed  using 
the  upper  and  lower  bound  theorems  of  plastic  analysis,  (20).  The 
lower  bound  theorem  states  that  if  a  yield  surface  is  derived  from  a 
safe  and  statically  admissible  stress  state  it  will  be  enclosed  by  or 
coincide  with  the  correct  yield  surface.  The  upper  bound  theorem 
states  that  if  the  states  of  stress  on  a  yield  surface  correspond  to 
assumed  mechanisms,  that  yield  surface  encloses  or  coincides  with  the 
correct  yield  surface.  These  theorems  may  be  combined  to  form  a 
uniqueness  theorem  which  states  that  a  yield  surface,  derived  from  a 
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safe  and  statically  admissible  stress  state,  and  corresponding  to  a 
kinematically  admissible  mechanism,  is  the  one  correct  yield  surface. 

The  yield  condition,  formulated  in  terms  of  normalized  stress 
resultants,  involves  only  those  stress  resultants  which  are  contained 
in  the  expression  for  the  strain  energy,  (16).  Referring  to  the  pre¬ 
vious  discussion  on  approximations  in  the  mixed  torsion  problem,  it  is 
assumed  that  the  effects  of  Saint  Venant  torsion  are  negligible.  The 
strain  energy  of  the  bar,  U  ,  of  Eq.  4,  reduces  to 


2U  *  NW'  +  M  U"  -  M  V"  +  M  d>" 
s  y  x  or 


(8) 


Hence,  the  yield  condition  may  be  expressed  in  terms  of  the  dimension¬ 
less  parameters 

(9a) 

(9b) 

(9c) 

(9d) 

in  which  N,  M  ,  M  ,  and  M  are  the  previously  defined  stress  resultants 
x  y  ui  r  7 

shown  in  Fig.  3  and  N  ,  M  ,  M  ,  M  are  their  fully  plastic  values. 

p  xp  yp 

Since  shear  stresses  have  been  neglected,  the  von  Mises  yield  condition 
reduces  to 

9  9 

(10) 


n  = 

=  N/N 

P 

m 

*  M  /M 

X 

x  xp 

m 

=  M  /M 

y 

y  yp 

m 

=  M  /M 

U) 

U)  cup 

a  2  n2 

a  =  a 
zz  y 


in  which  is  the  yield  stress  of  the  material  under  uniaxial  stress. 
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hypothesis  is  not  valid  for  the  cross  section  as  a  whole.  It  is  assumed, 
however,  that  Navier's  hypothesis  does  hold  for  each  element  of  the  cross 
section.  This  implies  that  there  is  one  neutral  axis  contained  in  each 
element  of  the  cross  section. 

The  equations  for  the  yield  surface  may  be  simplified  by  the  following 
approximations  for  thin  walled  sections  (19); 

1)  the  neutral  axis  passes  perpendicularly  through  any  element 
of  the  section,  as  shown  in  Fig.  5,  or, 

2)  if  the  neutral  axis  lies  wholly  within  an  element,  it  is 
parallel  to  the  long  dimension  of  the  element.  Fig.  5. 

To  illustrate  the  derivation  of  yield  surface  equations  for  a  thin 
walled  open  cross  section  including  warping  stresses,  consider  the 
Z-section  in  Fig.  6.  The  dimensions  of  the  section  which  has  an  enforced 
center  of  rotation  at  P  are: 

b  =  width  of  the  flanges; 

h  ■  depth  of  section,  center  to  center  of  the  flanges; 

A  =  area  of  one  flange;  A  =  area  of  web; 
t  w 

t  *  thickness  of  web  and  flanges;  c  *  A  /A,  =  h/b; 

w  f 

h  *  distance  from  P  to  center  of  top  flange  where  typically  h/h  '  1. 

The  fully  plastic  stress  resultants  are  found  to  be: 


■  Vt<2  +  c) 

(11a) 

W  '  *  °yAfh<4  +  '> 

(lib) 

[  *0A.b 

yp  y  f 

(11c) 

l  -  ^  0  A,  bh 
wp  2  y  f 

(lid) 

13 


In  general,  a  yield  surface  has  several  elements  each  related  to 
generic  neutral  axis  positions.  The  approximations  regarding  neutral 
axis  location  within  an  element  allow  a  reduction  in  the  number  of 
surface  elements  comprising  the  yield  surface  of  a  Z-section  to  four.  They 
are  shown  in  Fig.  7.  However,  since  the  section  is  not  symmetric,  four 
complementary  cases  are  necessary  for  complete  determination  of  the  yield 
surface.  In  these  complementary  cases,  indicated  in  Fig.  8,  the  generic 
neutral  axis  locations  of  Fig.  7  are  unchanged  but  the  sign  of  the 
yield  stress  is  reversed. 

The  stress  resultants  corresponding  to  the  stress  distri¬ 
bution  of  Fig.  9  are: 


N  *  oyAf(2  +  c  _  2(8  +  yc  +  a)) 

(12a) 

Mx  *  oyAfh(yc(y  -  1)  -  a  +  3) 

(12b) 

My  -  ayAfb(82  -  28  +  *2) 

(12c) 

M  *  2  °yAfbS(1  "  2a2) 

(12d) 

in  which  a,  8,  and  y  are  dimensionless  parameters  that  define  the 
neutral  axis  location  in  each  element  of  the  cross  section. 

The  corresponding  normalized  stress  resultants,  found  by  substi¬ 
tuting  Eqs.  12  and  11  into  Eqs.  9,  are: 

„ .  i  -  2tf,vc  *  a> 

2  +  c 

mx  “  4+^  (Yc(Y  "  l)  ~  a  +  S) 

2  2 
my  »  8  -  28  +  a 

m  -  1  -  2a2 

u> 


(13a) 

(13b) 

(13c) 

(13d) 


14 


Eliminating  the  parameters  a,  6,  and  y  from  Eqs.  13  ,  the  equation 
of  the  element  of  the  yield  surface  for  this  case  is  obtained  as 


Consider  the  stress  distribution  shown  in  Fig.  10.  In  this  case, 
a  neutral  axis  lies  within  the  bottom  flange.  The  values  of  n  and  m 

x 

are  insensitive  to  the  angle  8  since  it  is  necessarily  small.  The 
normalized  warping  moment,  m  ,  is  independent  of  0 ,  and  the  bending 
moment  about  the  y-axis,  m^,  is  an  unknown  function  of  0  and  hence 
indeterminate.  The  normalized  stress  resultants  are: 


m 


m 


m 


oj 


c  +  2(-l  +  B  +  a) 
2  +  c 

. 

4  +  c 


indeterminate 

2 

1  -  2a 


(15a) 

(15b) 

(15c) 

(15d) 


The  equation  for  this  element  of  the  yield  surface,  obtained  from  Eqs. 


15  ,  is  independent  of  m^  and  may  be  written 


(2  +  c)  n  +  ~  - 7  C  m  +  ^3(1  -  m  )*  +  c  -  2  *  0 
2  x  ’  w 


(16) 


Similarly,  equations  for  the  other  elements  were  found.  They  are 
shown  in  tabular  form  in  Figs.  7  and  8.  The  yield  surface  constructed 
from  the  various  surface  elements  is  illustrated  in  Fig.  11  for  the  Z- 
stiffener  of  Fig.  4.  A  two  dimensional  section  in  the  n  -  m^  plane 


is  shown  for  m 


m 


0.  The  range  of  applicability  of  each  surface 


element  is  indicated  on  the  figure.  Isometric  plots  of  the  vield 


DETERMINATION  OF  UPPER  BOUND-  -  YIELD  MECHANISMS  -  UNIQUENESS 

The  state  of  stress  in  the  bar  may  be  expressed  in  terms  of  the  stress 
resultants,  which  may  be  considered  components  of  a  matrix  denoted  by  { Q} : 

{Q}  =  {  N,  M  ,  M  ,  M  }  (17) 

y  x  oo 

The  corresponding  matrix  of  bar  strains,  {q} ,  is 

{q}  -  {  W',  0",  -  V",  4>"  }  (18) 


The  yield  surface  is  a  function  of  the  stress  resultants  such  that,  at 
yield, 


F  (».  My,  Mx,  «J  -  0 


(19) 


It  has  been  assumed  that  the  incremental  plastic  strain  is  directed 
outward,  normal  to  the  yield  surface.  This  associated  flow  rule  may 
be  stated  as 


.  ,  3F 


(20) 


in  which  A  is  an  arbitrary  positive  proportionality  factor  and  6q^ 
denotes  a  component  of  incremental  plastic  strain. 

From  Eq.  20  it  follows  that  the  strain  increments  corresponding 
to  stress  points  on  the  yield  surface  satisfy 


6q. 


3Q- 

3Qi 


(21) 


in  which  the  derivatives  on  the  right  hand  side  are  obtained  from 
Eq.  19  . 


To  illustrate  the  determination  of  such  strain  increment  rates, 
consider,  e.g.,che  lower  bound  stress  distribution  of  Fig.  10.  The 
yield  surface  equation  for  this  stress  state  is  given  by  Eq.  16. 


17 
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Comparison  of  Eqs.  22  and  24  shows  that  the  stress  state  on  this 
element  of  the  lower  bound  yield  surface  corresponds  to  a  yield  • 
mechanism.  Thus,  the  uniqueness  theorem  of  plastic  analysis  is 
satisfied  and  the  yield  surface  element,  given  by  Eq.  16,  is  the 
correct  one. 

In  a  similar  manner,  the  uniqueness  theorem  can  be  shown  to  be 
satisfied  for  the  remaining  surface  elements  of  Figs.  7  and  8. 

Hence,  the  complete  yield  surface,  defined  by  the  individual  surface 
elements,  is  the  exact  yield  surface.  The  yield  mechanisms,  derived 
from  kinematics,  for  each  surface  element  are  listed  in  Table  2. 
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ACCURACY  AND  IMPLEMENTATION  OF  THE  PROPOSED  FORMULATION 


The  accuracy  of  the  proposed  formulation  rests  on  two  approxi¬ 
mations.  The  first  involves  the  elimination  of  Saint  Venant  torsion 
from  the  yield  condition.  The  second  involves  assumptions  on  neutral 
axis  location  within  the  elements  of  the  cross  section. 

Referring  to  the  earlier  discussion  on  approximations  to  the 
mixed  torsion  problem,  it  is  evident  that  as  the  importance  of  Saint 
Venant  torsion  increases,  the  proposed  formulation  will  be  in¬ 
creasingly  in  error. 

The  assumption  that  a  neutral  axis  passes  perpendicularly 
through  an  element  of  the  cross  section  leads  to  artificial  ridges 
in  the  yield  surface.  These  ridges  lie  outside  the  correct  yield 
surface.  This  approximation  was  discussed  in  detail  in  (19),  and 
has  been  shown  to  be  negligible  in  elasto-plastic  analysis  of 
thin  walled  sections. 

The  technique  used  in  the  development  of  the  yield  surface  for 
a  Z-section  including  warping  stresses  can  be  applied  to  any  thin- 
walled  open  cross  section.  Other  common  types  used  as  ring-stiffeners 
for  cylindrical  shells  are  T  and  L  sections. 

The  proposed  yield  surface  equations  are  suitable  for  incor¬ 
poration  into  the  existing  nonlinear,  elasto-plastic  finite  element 
program,  EPSA  (2).  The  present  stiffener  formulation  in  the  code 
accounts  for  plastic  response  in  terms  of  only  two  stress  resultants, 
namely  the  axial  force  and  bending  moment  about  the  strong  axis  of  the 
stiffener.  The  proposed  formulation  will  extend  the  structural 
modeling  capability  for  the  class  of  problems  such  as  that  shown 
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in  Fig.  4,  and  lead  to  an  increase  in  the  accuracy  of  structural 
response  predictions  .especially  in  regions  of  large  plastic  strain 
and  deformation. 

CONCLUSION 

A  procedure  has  been  described  for  deriving  yield  surface 
equations  for  thin-walled  open  cross  sections  with  an  enforced  center 
of  rotation  for  cases  in  which  Saint  Venant  torsion  is  small  compared 
to  warping  torsion.  Equations  have  been  derived  for  the  specific 
case  of  a  Z-section.  These  equations  have  been  shown  to  satisfy  the 
uniqueness  theorem  of  plastic  analysis  thereby  furnishing  the 
correct  yield  surface  within  the  limits  of  validity  of  the  assumptions 
made.  The  applicability  and  accuracy  of  the  proposed  equations 


are  indicated. 
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NOTAT ION 


1 


The  following  svmbols  are  used 


in  this  paper: 


A.  =  cross  sectional  area  of  flange; 

b  =  width  of  flange; 

c  =  ratio  of  area  of  web  to  area  of  flange; 

d  =  characteristic  length; 

E  =  Young's  modulus; 

F  =  function  defining  yield  surface 

G  =  shear  modulus 

depth  of  section  center  to  center  of  flanges; 
distance  from  pole  to  center  of  top  flange; 

I  =  warping  constant; 

J  =  torsion  constant; 

L  =  length  of  bar; 

M  =  bending  moment  about  x-axis; 

m  =  normalized  bending  moment  about  x-axis; 

x 

M  “  fully  plastic  bending  moment  about  x-axis; 

M  =  bending  moment  about  y-axis; 

V 

m^  =  normalized  bending  moment  about  y-axis; 

M  =  fully  plastic  bending  moment  about  y-axis; 

M  =  warping  moment; 


m  = 

JU 

M  = 

•jJ  p 

N  = 
n  = 

N  = 

_P 

n-s-z  = 
P  = 

Q  = 
q  = 

r  = 

T  = 
s 

t  * 

u.v.w,- 

u  * 

s 

v  * 


x,y,z  » 

a.S.Y  * 


5qi 


normalized  warping  moment; 

fully  plastic  warping  moment; 

axial  force; 
normalized  axial  force; 
fully  plastic  axial  force; 

coordinate  system  defined  at  middle  surface  of  cross  section ; 
pole  or  center  of  rotation; 
matrix  of  stress  resultants; 
matrix  of  oar  strains; 

perpendicular  distance  from  pole  to  point  on  middle 
surface  of  section; 

Saint  Venant  torque; 

thickness  of  wall  of  cross  section 

rigid  body  displacement  of  pole  in  the  x,y,z,  directions; 
strain  energy  per  unit  length  of  bar; 

axial  displacement  of  point  on  the  middle  surface  of 
the  cross  section 
cartesian  coordinates; 

dimensionless  parameters  locating  neutral  axis  in  cross 
section; 

increment  of  plastic  bar  strain; 
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TABLE  1  -  LOCATION  OF  NEUTRAL  AXES  FOR  SURFACE  ELEMENT  2 
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TABLE  2  -  YIELD  MECHANISMS 
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FIG-1  CROSS  SECTION  DISPLACEMENTS 


FIG.  2  STRESS  COMPONENTS 


FIG.  3  STRESS  RESULTANTS 
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FIG.  5  ASSUMED  NORMAL  STRESS  DISTRIBUTIONS 
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DEPARTMENT  OF  THE  NAVY  (Continued) 

Naval  Surface  Weapons  Center 
ATTN:  W.  Wishard 

ATTN:  Tech  Library  &  Info  Svcs  Br 

Naval  War  Col  lege 

ATTN:  Code  E-ll 

Naval  Weapons  Center 
ATTN:  Code  233 
ATTN:  Code  266,  C.  Austin 
ATTN:  Code  3263,  J.  Bowen 

Naval  Weapons  Evaluation  Facility 
ATTN:  G.  Binns 
ATTN:  Code  10 
ATTN:  Code  210 
ATTN:  R.  Hughes 

Naval  Weapons  Support  Center 

ATTN:  Code  70553,  D.  Moore 

New  London  Laboratory 

ATTN:  Code  4494,  J.  Patel 
ATTN:  Code  4492,  J.  Kalinowski 

Newport  Laboratory 
ATTN :  Code  EM 

ATTN:  Code  363,  P.  Paranzino 


Naval 

Ocean  Systems  Center 

Ofc  of  the  Deputy  Chief  of  Naval  Ops 

ATTN 

Code  013,  E. 

Cooper 

ATTN 

OP  987 

ATTN 

Code  4471 

ATTN 

NOP  982,  Tac  Air  Srf  6  Ewdev  Div 

ATTN 

NOP  981 

Naval 

Postgraduate  School 

ATTN 

NOP  654,  Strat  Eval  5  Anal  Br 

ATTN 

Code  69NE 

ATTN 

OP  Q98T8 

ATTN 

Code  1424,  Library 

ATTN 

OP  982E,  M.  Lenzini 

ATTN 

Code  69SG,  Y. 

Shin 

ATTN 

OP  957E 

ATTN 

NOP  953,  Tac  Readiness  Div 

Naval 

Research  Laboratory 

ATTN 

OP  37 

ATTN 

Code  3403,  R. 

Bel  sham 

ATTN 

OP  225 

ATTN 

Code  8440,  G. 

O'Hara 

ATTN 

OP  03EG 

ATTN 

Code  6380 

ATTN 

OP  21 

ATTN 

Code  8100 

ATTN 

NOP  951,  ASW  Div 

ATTN 

Code  8301 

ATTN 

OP  60505 

ATTN 

Code  8406 

ATTN 

OP  981N1 

ATTN 

Code  2627 

ATTN 

OP  223 

ATTN 

Code  8445 

ATTN 

Code  8404,  H. 

Pusey 

Office  of  Naval  Research 

ATTN 

Code  474,  N.  Perrone 

Naval 

Sea  Systems  Command 

ATTN 

SEA-033 

Strategic  Systems  Project  Office 

ATTN 

SEA- 323 

ATTN 

NSP-272 

ATTN 

SEA-06J,  R.  Lane 

ATTN 

NSP-43 

ATTN 

SEA-09G53 

ATTN 

NSP-273 

ATTN 

SEA- 55X1 

ATTN 

SEA-08 

DEPARTMENT  OF  THE  AIR  FORCE 

ATTN 

SEA-0351 

ATTN 

SEA-9931G 

Air  Force 

nstitute  of  Technology 

ATTN :  Coitmander 


Naval  Surface  Weapons  Center  ATTN:  Library 


ATTN 

Code  F34 

ATTN 

Coe-  R13 

Air  Force  Systems  Command 

ATTN 

Code  :'10 

ATTN:  XW 

ATTN 

Code  U401,  M.  Klelnerman 

ATTN 

Code  R14 

Assistant  Chief  of  Staff 

ATTN 

Code  F31 

Intel  1 igence 

ATTN 

Code  R15 

ATTN:  IN 
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DEPARTMENT  Of  THE  AIR  FORCE  (Continued) 


DEPARTMENT  OF  ENERGY  CONTRACTORS  ( Conti nued) 


Air  Force  Weapons  Laboratory 

L.os  Alamos 

National  Laboratory 

ATTN 

NTES-G,  S.  Melzer 

ATTN 

R.  Whitaker 

ATTN 

NTE,  M,  Plamondon 

ATTN 

MS530,  G.  Spillman 

ATTN 

NTES-C.  R.  Henny 

ATTN 

Reports  Library 

ATTN 

SUL 

ATTN 

MS634,  T.  Dowler 

ATTN 

NTED 

ATTN 

R.  Sanford 

ATTN 

MS670,  J.  Hopkins 

Ballistic  Missile  Office 

ATTN 

DEB 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Deputy  Chief  of  Staff 

Applied  Research  Associates,  Inc 

Research,  Development,  &  Acq 

ATTN 

D.  Piepenburg 

ATTN 

AFROqi 

ATTN 

R.  Steere 

Applied  Research  Associates.  Inc 

ATTN 

B .  Frank 

Deputy  Chief  of  Staff 

Logistics  A  Engineering 

BOM  Corp 

ATTN 

LEEE 

ATTN 

T.  Neighbors 

ATTN 

A.  Lavagnino 

Foreign  Technology  Division 

ATTN 

Corporate  Library 

ATTN 

NIIS  Library 

ATTN 

TQTO 

California 

Institute  of  Technology 

ATTN 

SDBG 

ATTN 

T.  Ahrens 

ATTN 

SDBF,  S.  Spring 

Cal  ifrmia 

Research  A  Technology,  Inc 

Rome  Air  Development  Center 

ATTN 

M.  Rosenblatt 

ATTN 

RBES,  R.  Mair 

ATTN 

S.  Schuster 

ATTN 

Commander 

ATTN 

Library 

ATTN 

TSLD 

ATTN 

K.  kreyenhagen 

Stragetic  Air  Command 

Columbia  University 

ATTN 

NRI-STINFO  Library 

ATTN 

H.  Bleich 

ATTN 

F.  Dimaggio 

DEPARTMENT  OF  ENERGY 

University 

of  Denver 

Albuquerque  Operations  Office 

ATTN 

Sec  Officer  for  J,  Wisot' 

ATTN 

CTID 

Electric  Power  Research  Institute 

Department  of  Energy 

ATTN 

G.  Sliter 

ATTN 

OMA/RDAT 

Electro-Mech  Systems,  Inc 

Nevada  Operations  Office 

ATTN 

R.  Shunk 

ATTN:  Doc  Con  for  Technical  Library 

General  Dynamics  Corp 


OTHER  GOVERNMENT  AGENCIES 

ATTN 

J.  Miller 

ATTN 

J.  Mador 

Central  Intelligence  Agency 

ATTN 

M.  Pakstys 

ATTN: 

OSWR/NED 

ATTN: 

OSR/SE/F 

Kaman 

AviDyne 

ATTN 

R.  Ruetenik 

NASA 

ATTN 

G.  Zartarian 

ATTN: 

F.  Nichols 

ATTN 

Library 

ATTN: 

R.  Jackson 

ATTN 

N.  Hobbs 

NATO 

kaman 

Sciences  Corp 

ATTN 

Library 

NATO  School 
SHAPE 

ATTN 

F.  Shelton 

ATTN: 

US  Oocuments  Officer 

Hainan 

Sciences  Corp 

ATTN:  D.  Sachs 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

Kaman 

Tempo 

University  of  California 

ATTN 

DAS I AC 

Lawrence  Livermore  National  Lab 
ATTN:  S.  Erickson 
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OEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Karagozian  and  Case 

ATTN:  J.  Karagozian 

Lockheed  Missiles  A  Space  Co,  Inc 

ATTN:  Technical  Information  Center 
ATTN:  T.  Geers 
ATTN:  B.  Almroth 

Lockheed  Missiles  A  Space  Co,  Inc 
ATTN:  TIC-Library 

MAT  Company 

ATTN:  0.  McNaight 

Management  Science  Associates 
ATTN:  K.  Kaplan 

McDonnell  Douglas  Corp 
ATTN:  R.  Halprin 

NKF  Engineering  Associates,  Inc 
ATTN:  R.  Belsheim 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 

Pacifica  Technology 
ATTN:  R.  Bjork 
ATTN:  G.  Kent 
ATTN:  A.  Kushner 

Physics  Applications,  Inc 
ATTN:  C.  Vincent 

Physics  International  Co 


ATTN: 

L. 

Behrmann 

ATTN: 

F. 

Sauer 

ATTN: 

J. 

Thomsen 

ATTN: 

E. 

Moore 

ATTN: 

Technical  Library 

Science  Appl icatlons,  Inc 

ATTN:  Technical  Library 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued 

S-CUBED 

ATTN 

T.  Cherry 

ATTN 

R.  Sedgewick 

ATTN 

D.  Grine 

ATTN 

T .  Ri  ney 

ATTN 

Library 

ATTN 

K,  Pyatt 

ATTN 

T.  McKinley 

SRI  International 

ATTN 

G.  Abrahamson 

ATTN 

W.  Wilkinson 

ATTN 

A.  Florence 

Teledyne  Brown  Engineering 

ATTN 

J.  Ravenscraft 

Tetra  Tech,  Inc 

ATTN 

L.  Hwang 

TRW  Electronics  4  Defense  Sector 

ATTN 

P.  Bhuta 

ATTN 

A.  Feldman 

ATTN 

N.  Lipner 

ATTN 

Technical  Information  Center 

ATTN 

D.  Jortner 

ATTN 

B.  Sussholtz 

TRW  Electronics  4  Defense  Sector 

ATTN 

P.  Dai 

ATTN 

F.  Pieper 

ATTN 

E.  Wong 

ATTN 

G.  Hu! cher 

Weidlinger  Assoc,  Consulting  Engrg 

ATTN 

J.  McCormick 

ATTN 

M.  Baron 

4  cy  ATTN 

R.  Daddazio 

4  cy  ATTN 

M.  Bleniek 

4  cy  ATTN 

F.  Dimaggio 

Weidlinger  Associates 

ATTN 

J.  Isenberg 

Southwest  Research  Institute 
ATTN:  A.  Wenzel 
ATTN:  W.  Baker 


AO 


